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Abstract: Effective solutions should be devised to handle the effects of shadow zones in 
Underwater Wireless Sensor Networks (UWSNs). An adaptive topology reorganization 
scheme that maintains connectivity in multi-hop UWSNs affected by shadow zones has 
been developed in the context of two Spanish-funded research projects. A mathematical 
model has been proposed to find the optimal location for sensors with two objectives: the 
minimization of the transmission loss and the maintenance of network connectivity. The 
theoretical analysis and the numerical evaluations reveal that our scheme reduces the 
transmission loss under all propagation phenomena scenarios for all water depths in 
UWSNs and improves the signal-to-noise ratio. 
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1. Introduction 

Underwater Wireless Communication Networks (UWCNs) are formed by sensors and Autonomous 
Underwater Vehicles (AUVs) interacting together to perform specific underwater applications such as 
collaborative monitoring or surveillance [1]. Communication quality in Underwater Wireless Sensor 
Networks (UWSNs) is very challenging due to the harsh characteristics of the underwater channel, 
such as high and variable propagation delays, limited bandwidth, high bit error rates, multipath 
phenomena and multipath fading [1]. In the extreme case, the spatially- variant underwater channel can 
cause the formation of shadow zones, which are time-variant areas where there is little signal 
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propagation energy due to the refraction of signals by the sound speed fluctuation [2]. When the sound 
speed has a negative gradient just beneath the surface [3], a shadow zone is formed because the 
acoustic rays are refracted downward. Refraction produces shadow zones that sound waves do not 
penetrate because of their curvature. The sea bottom can produce a shadow zone as well, when the 
rays are refracted upward. Shadow zones [3] can also appear beneath the mixed layer for a source 
located near the ocean surface because the acoustic energy is trapped in the surface duct (see Figure 1). 
The shadow zone is usually bounded by the lower boundary of the surface duct and the limiting ray. 
Shadow zones can also appear between convergence zones. If the source is located at the same depth 
of the underwater sound channel axis, the shadow zone will disappear [4]. We distinguish between 
shallow (depth up to 100 m) and deep water. In shallow water (order of 100 meters depth) and at 
ranges of 3 kms shadow zones appear [2]. In deep water (order of 1,000s of meters depth) and at 
ranges of 10s of kms shadow zones are formed [2]. 



Figure 1. Shadow zone formation beneath the mixed layer when sound velocity 
monotically decreases with depth. 

0 c 0 c c 




Shadow zones cause high bit error rates, losses of connectivity and dramatically impact 
communications performance. Some experiments show that for high frequencies signal levels are 
typically at least 40 dB less than those at the edges of the shadow zone [5]. For low frequencies the 
signal loss in the shadow zone is less severe; however, variations in received signal-to-noise ratios 
(SNRs) by as much as 10 dB have been observed on time scales of several hours [6]. Therefore, many 
authors recommend devising solutions to handle their effects [7,8]. 

For these reasons, in this paper we focus on communication reliability in the presence of shadow 
zones. We propose a distributed adaptive topology reorganization scheme that alleviates the effects of 
energy limitations and is able to maintain connectivity between sensor nodes in multi-hop 
three-dimensional UWSNs in the presence of shadow zones. Besides, it is able to estimate when the 
shadow zones have disappeared using double sensor units to restablish communication very quickly 
through the original acoustic wireless links. We study the effects of the proposed scheme in shallow 
(depth up to 100 m) and deep water. A two-path Rayleigh fading channel model and different 
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propagation phenomena are considered: shallow water, deep water with convergence zones, deep 
water with deep sound channel and shallow or deep water with shadow zones. 

According to the number of links in the three-dimensional UWSN affected by the shadow zone, 
three different cases have been introduced. In this paper we have extended our work in [9], where only 
one case was discussed. A mathematical model has been developed for each case to find the optimal 
placement for the sensor nodes whose ongoing communications are being disturbed. The two major 
objectives of the proposed mathematical model are: the minimization of the transmission loss and the 
maintenance of network connectivity. To the best of our knowledge this is the first paper that 
finds the optimal location for underwater sensor nodes affected by a shadow zone under different 
propagation phenomena. 

The theoretical analysis and the numerical evaluations reveal that the average transmission loss 
values are reduced significantly under all propagation phenomena for all water depths in UWSNs 
when the optimal locations are computed and the communication between sensor nodes is again 
restablished outside the shadow zone. The average SNR values have also improved significantly and 
are maintained for all frequencies. 

The remainder of the paper is organized as follows: In Section 2, we discuss the related work. In 
Section 3, we analyze our system model. In Section 4, we state the location optimization problems and 
propose nonlinear programming (NLP) formulations. In Section 5 we present our numerical results. 
Finally, we draw the conclusions in Section 6. 

2. Related Work 

Shadow zones represent a serious obstacle for good communication in UWSNs, because they cause 
network partition. Therefore, the proposal of solutions to handle their effects has been encouraged [7,8]. 
A distributed reactive shadow zone and delay aware (SZODAR) routing protocol for 
three-dimensional UWSNs was introduced in [10]. In SZODAR, the depth of the acoustic transceivers 
of sensor nodes is changed to avoid shadow zones, while the sensing modules are kept unchanged. 
During the connection establishment process, each node should send HELLO messages until it 
receives a routeReply or until the maximum number of allowed HELLO messages maxTrial is 
reached. If the number of HELLO messages is over maxTrial, the node is considered to be out of range 
or located in a shadow zone. In this case the node should increase its transmission power to the next 
power level stepwise up to the maximum and try to establish connection again at each step. Since 
UWSNs have very limited resources, the proposed connection establishment process consumes 
excessive power and can lead easily to energy depletion. In addition, it is not possible to determine if 
the node unable to establish connection using maximum power level is located in a shadow zone or is 
simply a node out of the range of the other nodes. Our proposed scheme has been designed to detect 
shadow zones based on the estimation of the transmission loss values, which is a not such an excessive 
energy consuming operation. 

In addition, in [10] the acoustic transceiver of the node out of range or located in a shadow zone is 
moved upwards or downwards according to a random probability, whereas the node periodically tries 
to establish connection with its neighbours. If the boundary for the operational depth is reached and the 
connection establishment process has failed, the node moves in the opposite direction. With this 
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arbitrary movement criterion this shadow zone awareness routing protocol consumes excessive power 
and incurs high delays and routing overhead. Moreover, the effects of different underwater 
propagation phenomena on the performance of the proposed protocol have not been analyzed. 

In our proposed scheme sensor nodes are double units operating as a single sensor; they are 
decoupled into two sensor nodes in the presence of a shadow zone. One sensor node remains in the 
same position and estimates, as opposed to SZODAR [10], when the shadow zone has disappeared. An 
optimal placement imposing a restriction in energy consumption is computed for the other sensor 
node; this node is moved directly to this optimal location away from the shadow zone. In this way, 
connectivity between both sensors inside the shadow zone is established by a wire, connectivity 
between neighbour nodes is maintained and latency as well as power consumption are reduced. The 
effects of the different propagation phenomena on the performace of the proposed scheme have been 
analyzed. Furthermore, connection can again be reestablished very quickly using the original wireless 
acoustic links if the shadow zones disappear. 

Several papers have proposed options related to location optimization of sensor nodes underwater. 
In [11], a placement scheme has been introduced to find the optimal placement of data collectors in 
UWSNs. In [12], a distributed node deployment technique for UWSNs to improve coverage and 
provide connectivity with the surface stations has been proposed and compared in terms of coverage 
and connectivity with the optimal centralized solutions from the literature. In [13], deployment 
strategies for two-dimensional and three-dimensional communication architectures for UWSNs have 
been developed to determine the minimum number of sensors to be deployed to achieve optimal 
sensing and communication coverage according to the application requirements. However, these 
contributions don't analyze which is the optimal location for underwater sensor nodes affected by a 
shadow zone under different propagation phenomena. 

3. System Model Analysis 

We consider a three-dimensional underwater sensor network for environmental monitoring 
(see Figure 2). In this architecture, a certain number of underwater sensors are deployed to cover a 
large area of a lake or ocean. Metal pieces at the seabed and thin ropes are used to attach the sensor 
nodes and distribute them at different depth levels [1]; they are deployed at different water depths to 
observe water quality parameters since water characteristics vary mostly in vertical direction. An 
electronically controlled engine at the sensors is used to change their depths [1]. This deployment 
strategy requires building a tree-like hierarchical multi-hop routing topology. The sensed data is sent 
to the surface station, which floats on the water surface and uses long distance radio communication to 
send this data to an onshore station. The onshore station is located at the monitoring center and 
performs further analysis of the data. 
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Figure 2. Architecture for a 3D underwater sensor network. 




4. Problem Formulations and Optimization Models 

We consider the example of the tree-like topology shown in Figure 2, where sensor nodes sense 
data and sent it upwards to the surface station using multi-hop routing. We assume that the distributed 
sensor nodes are double units. The cost of implementing double sensors will be limited because 
software modems will be coupled with generic microphones and speakers that are built-in to create 
off-the-shelf sensor modules [14-16]. In this way, the need for specialized communication hardware 
will be eliminated and the system cost will be reduced. The cost per node ($) for hardware modems 
that use especialized hardware is around 10 4 $, whereas the cost of software modems that use generic 
hardware (off-the-shelf configuration) is around 10 2 $; this cost decrease facilitates the deployment of 
sensor nodes to form underwater acoustic sensor networks. The double sensor nodes operate as a 
single sensor when the underwater communication is reliable. This fact changes in the presence of 
shadow zones. 

Each intermediate node Si should be able to detect if the quality of the uplink towards one-hop 
neighbor s i+ \ is affected by the presence of a shadow zone computing the transmission loss TLu+i [3]. 
If TLi i+ i(t + A^) = TLi i+ i(t) + TL shadow zone and this value persists during a period longer than T to 
eliminate insignificant variations, the sensor node assumes that the uplink is located in a shadow zone 



Sensors 2009, 9 



8689 



because TL t i+ \(t + At) has been increased by the transmission loss value TL shadow zone inside the 
shadow zone. 

We distinguish between three different cases according to the number of links in the underwater 
sensor network affected by the shadow zone: 

• Case 1 : Only the uplink i+ \ between s t and s i+ \ is located in the shadow zone. 

• Case 2: The uplink lu+\ between Si and s i+ \ and the uplinks l(i-\)ki between S(i.\)k 9 Vk e K and 
Si are located in the shadow zone. K represents the set of neighbours connected to Si through 
uplinks using hierarchical routing. 

• Case 3: The uplink between s t and s i+ \ and the uplinks l(i-\)ki between S(i-\)h Vk e K and 
Si, l(i-2)k i-\ between s^h Vk e K and a particular s (i .i)h /(*-*)* i-x+i between s (/ _ x )£, 
\/k g K, Vx e and a particular S(i. x +\)k are located in the shadow zone. 

Each case will be solved differently. In all cases we assume that the sensor node s t knows the 
location and dimensions of the shadow zone. The sound speed profile is affected by the time of the 
day, weather, depth, variations in the water temperature and salinity. With the aid of a forecast 
model [8] and the estimation of the sound speed profile, the different underwater propagation 
phenomena can be predicted, including the presence of shadow zones. The location of shadow zones 
can be performed studying the trajectory of the acoustic rays under the different underwater 
propagation phenomena [3]. Autonomous Underwater Vehicles (AUVs) can analyze and inform 
underwater sensors periodically about the existence and location of shadow zones. 

We assume that sensor nodes are located in a grid (see Figure 6). The axes of the coordinate system 
are set at the water surface. Sensors know their fixed two coordinates (x,y). They are also aware of 
their depth (z coordinate), which varies with time depending on the networks tasks. The z coordinate is 
zero at surface and increases with depth. Sensor nodes also know the location of their neighbours 
because they exchange their depth information. 

Case 1: Only the uplink la+\ between Si and Si+\ is located in the shadow zone. 

In this case under the presence of a shadow zone the sensor node s t is uncoupled into two sensor 
nodes s t ->\ and Si-& as shown in Figure 3. These sensor nodes will be connected to each other inside 
the shadow zone by a wire to maintain robust communication. The node s^\ remains in the same 
location as the node Si was and continues sensing data at the required depth level. We should now 
determine the optimal location for the node Si-& We formulate this problem as a Nonlinear 
Program (NLP). 

We introduce the following notation: 

• Ui+\ = (x i+ i, yt+h Zi+i) is the location of node s i+ \. 

• u^2 = (xu yu Zi^i) is the new location of node s^2. 

• His the sea depth. 

• Rmax is the maximum transmission range. 

• TL T h is the detection threshold. 

• 7i->2 i+i refers to the transmission range of the acoustic link between nodes 5,-_>2 and s i+ \ 
expressed in meters. 
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a represents the absorption coefficient and has the units dB/Km. 

X is a variable used to compute the transmission loss with the value % = 10 for shallow water 
(cylindrical spreading) and %= 20 for deep water (spherical spreading) [3]. 
A is called the transmission loss anomaly and summarizes the combined effects of several 
complex factors such as multiple path propagation due to the variations of the speed with 
temperature, depth and salinity, refraction effects, diffraction and scattering of sound by 
particles, bubbles and plankton within the water column. Detailed formulas of the multipath 
propagation loss according to the different propagation phenomena (surface reflections, 
bottom reflections, convergence zones, deep sound channel, etc.) can be found in [3]. It is 
given in dB. 

Wsz-start = yu Zszstari) and w sz -end = yu z sz - end ) are the coordinates of the points 
delimiting the shadow zone at (x h y t ) with z sz - start < z sz - end and height z sz - end - z sz - start . 
E Bw is the energy consumed moving the sensor s^2> 
E T h is the energy consumption threshold. 



Figure 3. Scenario in the 3D UWSN with a shadow zone (Case 1). 
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The problem can be formulated as follows: 



PI: Location Optimization Problem 
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Min : TL M M = z log r^ 2 M + a.r^ 2 i+l . 1 0" 3 + ^ 
Subject to : 

^ 2i+1 =|K + i-^ 2 | (i) 

r,_> 2 <+1 < R MAX (2) 
TL^ 2M <TL Th (3) 

Z z^2 < Z sz-start (^) 
< ^7% ( 5 ) 

The objective function of problem PI aims at finding the best location for s^2 that minimizes the 
transmission loss and preserves the connectivity with node s i+ \ preventing from network partition. 

The absorption coefficient a is computed as derived in [17,18] as a function of the frequency /, the 
salinity S, the temperature r, the water ph pH and the depth z. We consider a two-path Rayleigh fading 
channel model, where each path is Rayleigh distributed. More specifically, the envelope of the signal 
from each path in the two-path model is modeled as an independent Rayleigh distributed random 
variable, a u i e {1,2} [19]. Consequently, for each path, the received energy per bit per noise power 
spectral density is given by, y/= a 2 .(Eb/No) which has a distribution, fT(y/)=(l/y/o).exp(-y//y/o) where 
y/o = E[a 2 ].(Eb/No) and Ej/No can be found from the SNR of the channel. Since the 2-path Rayleigh 
model does not have a closed- form expression for SNR distribution, it is found through simulations [19]. 
The value of a is used to find the transmission loss and the SNR of the channel. 

Constraint (2) imposes that the transmission range should be lower than a threshold to ensure 
connectivity. Constraint (3) expresses that the transmission loss should be lower than a transmission 
loss threshold, that is, the maximum propagation loss for properly receiving the transmitted signal. 

The SNR of an emitted underwater signal at the receiver can be expressed in dB by the passive 
sonar equation [3]: 

¥ , 4B =SL-TL-NL + DI> y/ Th (6) 

where y/Th has been defined as the detection threshold, SL is the source level, TL is the transmission 
loss, NL is the noise level and DI is the directivity index. 

The signal level SL is related to the intensity I t and hence to the transmission power of the 
transceiver for shallow water as follows: 



I t = — ^ 9 sz = 101og 

2n \m z 



0.67xl(T 18 y 



(7) 



where P t is the transmission power in watts and z is the depth in meters 
In deep water the equation (7) becomes: 

P ( 
I t = , SZ = 101og 

4tt Im 



0.67 xl(T 18 



(8) 



In our analysis we assume a fixed detection threshold for the received SNR. Based on this 
assumption we define the transmission loss threshold. It is given by [20]: 
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NL-y/ Th 

where SL is fixed and NL is also fixed for a particular operating frequency. 
We consider five different channel models: 

• Shallow water. 

• Deep water with convergence zones. 

• Deep water with deep sound channel. 

• Shallow water with shadow zones. 

• Deep water with shadow zones. 
Details of these models can be found in [3]. 

Constraint (4) imposes that the sensor node 5,-_>2 should be located outside the shadow zone. 
Constraint (5) states that the energy consumed by the electronically controlled engine to move the 
sensor to the optimal position should be lower than a threshold. It is defined as [21]: 

E Bw = P Bw .T Bw + A.v s .T Bw (10) 

where P Bw is the consumed power to move the sensor upwards or downwards and depends on the 
power required to drive the electronics, T Bw is the time required for the movement, v s is the moving 
speed and X is a constant depending on factors such as the drag force, sensor size and shape, friction of 
the motor, etc. With the proposed mechanism some extra energy will be consumed to locate the 
shadow zone and to relocate the sensor to the required position, but communication in the shadow 
zone will be possible. 

Case 2: The uplink 1,- i+ \ between Si and s i+ \ and the uplinks l{j.\)k / between •$(,•_!)& V£ e K and Si are 
located in the shadow zone. K represents the set of neighbours connected to s t through uplinks using 
hierarchical routing. 

In this case under the presence of a shadow zone the sensor node Si is uncoupled into two sensor 
nodes s^\ and 5,-_>2, as shown in Figure 4. These sensor nodes will be connected to each other inside 
the shadow zone by a wire to maintain robust communication. Now we should determine the optimal 
location for the nodes s^i and 5,-_>2 using NLP. Since the topology is tree-like, we assume several 
nodes S(i.\)k 9 Vk e K forward their data messages towards node s t . Our objective consists of minimizing 
the transmission loss while preventing from network partition. For this purpose the connectivity 
between nodes 5,-_>i and S(i.\)h e K on the one hand, and the connectivity between nodes 5,-_>2 and 
Si+\ on the other hand should be maintained. The optimal location for the node s^2 is obtained solving 
the PI optimization problem already introduced in Case 1. We should now obtain the optimal 
placement for the node s^i. 

We define: 

• um = (xu y u z,-_>i) is the new location of node s^\. 

• U(i-i)k= (x(i-i)h y{i-\)k, Z(i-i)k), V£ g K is the location of node s (i .i)k • 
Given : x. , x {i _ l)k , y t , y {i _ l)k , z {i _ A , H , R MAX TL Th , 

a -> %•> A? Z sz-start ' Z sz-end ' ^ Bw ' ^Th 
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Find: z. x g[0 9 H] 9 n 



i? 4 



Min max, 7Z( M) ^ 



Subject to : 



r (M), ^1 



ykGK 



r {i _ x)k ^<R MAx ykeK 



TL (i _ l)ki ^<TL Th ykeK 
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(11) 
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(13) 

(14) 
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Figure 4. Scenario in the 3D UWSN with a shadow zone (Case 2). 
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P2: Location Optimization Problem 

The objective function of problem P2 aims at finding the best location for s^\ that minimizes the 
transmission loss of the link ,■ with the highest transmission loss value and preserves the 
connectivity with nodes sy-iyt preventing from network partition. 

Constraint (12) imposes that the transmission range should be lower than a threshold to ensure 
connectivity. Constraint (13) expresses that the transmission loss should be lower than a transmission 
loss threshold. Constraint (14) imposes that the sensor node should be located outside the shadow 
zone. Constraint (15) states that the energy consumed by the electronically controlled engine to move 
the sensor to the optimal position should be lower than a threshold. 
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Case 3: The uplink la+\ between Si and Si+\ and the uplinks l{j.\)ki between S(i-\)k, V£ e K and su l(i-2)ki-i 
between s^k, Vk e K and a particular s^.\yh l(i- x )ki-x+\ between S(i- X )k 9 V£ e K, Vx e N and a 
particular S(i. x +\)k are located in the shadow zone. 

In this case the sensor node located at the lowest depth in the shadow zone s t is uncoupled into two 
sensor nodes s^\ and s^ as shown in Figure 5. 

Figure 5. Scenario in the 3D UWSN with a shadow zone (Case 3). 
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The optimal placement for node s z -_>2 is obtained solving the PI optimization problem already 
introduced in Case 1. The node s^\ should move to (jc,-, yu z sz -end + where s g R + is a very small 
value, that is, s^\ moves vertically downwards outside the shadow zone. In this case, the connectivity 
with the nodes s^k, \/k e K through the uplinks can only be maintained if they also move outside the 
shadow zone. Therefore, the sensor nodes sy-iyk, \/k e K are uncoupled into two sensor nodes SQ^yk-^i 
and S(i-iyk-+2- The node S( z --i)^2 remains in the same location as node sp-i)*, \/k e K was and continues 
sensing data at the required depth level. We should obtain the optimal placement for the node S(i-iyt-+i- 
The data sensed by £(m)£_>2 is sent through a wire between S(m)£_>2 and S(i.\)k-^i to maintain 
communication inside the shadow zone. Generally speaking, the same process should be repeated for 
all the nodes down in the hierarchy s^ x yk 9 \/k e K, \/x e N located in the shadow zone. They are 
uncoupled into two sensor nodes S(i- x yc-+i and S(i- x yc-+2> The node SQ- x )k-+2 remains in the same location as 
node S(i- X )k, V£ e K, Vx g N was and continues sensing data at the required depth level. We should 
obtain the optimal placement for the node S(i- x yk-+i. The data sensed by s^ x yk-^2 is sent through a wire 
between S(i- x yc-+2 and s^xyt^i to maintain communication inside the shadow zone. 

If the node S(i. x +\)k-^i has already moved to a new location outside of the shadow zone and node 
does not know its depth z^+i^i, node s^. x ) k ^i can move to (x {i . x ) h y {i . x)h z sz _ end + s) and once 
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there establish connection with node so- x +i)^i to find out the depth and afterwards compute its optimal 
location using this information. 

Finally, the data sent towards the suface sink travels through the following nodes in the path using 
hierarchical routing: ^-^2, s (i . x)k ^u . s (i .i)k^u ^1, ^2, s i9 s i+ i,. . ., and the surface sink. 

Generally speaking, the optimal location of the node 1 can be determined as follows. We define: 

• U(i- X ) k -+i = (x(i. X )hy(i-x)hZ(i-x)k^i)yk e K, Vx e Nis the new location of node s (i . x)k ^i. 

• U(i- x+ \)k-+\ = (x(i- x n)k,y(i-xn)k,Z(i-x+i)k^i) is the new already established location of node sy. 
jc+i)it->i for a particular k e K outside the shadow zone. We consider that node S(i- x +i)k was 
located originally in the shadow zone. 

• U(i- X .i)k= (x(i. X A)k,y(i-xA)k,Z(i-xA)k)yk g K f Vx g N is the location of node S(i. x .i)h We consider 
that node S(i. x -i)k is not located in the shadow zone. 

P3: Location Optimization Problem 

Given : x {i _ A , x {i _ x+A , x (i _ x _ l)k , y (i _ x)k , y (i _ x+l)k , y (i _ x _ l)k 

Z (i-x+l) k ^l ? Z (i-x-l) k ? H> ^MAxJ^Th ? ^? >3f> ^? 

Z sz-start ? Z sz-end ? ? ^77* 

Find : zp.^, e [0, //], r^^, (i . x+i)l ^ e 2T 

r (i-x-l) t (i-*) t -H G ^ 

Objective 1 : 

Mfo: ^o-x)^i o-x + i)^i 
Or Objective 2 : 

Min max, (2^4 (i . x)k ^,TL (i _ x)k ^ (i . x+l)i ^) 
Subject to : 



r ('-x) k ^l («-*+l)t-»l 



r (i-x-l) t (/-*)* 



,VkeK,VxeN (16) 
VkeK,VxeN (17) 



(! ._ x+1)t ^ < /? MAY , V* e ^,Vx g TV (18) 

< iW V£ eK,VxeN (19) 

?V*) t -i < 7Z ra ,V£ e tf, Vx e N (20) 

7V*_ l)t ('-*).-»! < TL ™> yk zK,\/xeN (21) 

Z (i-x) k -+l > Z sz-end 

(23) 

We consider two different objectives. The objective function 1 of problem P3 aims at finding the 
best location for S{i- X )h-+\ that minimizes the transmission loss of the link l(i- x )k-^i(i-x+i)k^\ and preserves 
the connectivity with node S(i. x +i)k^i preventing from network partition. The objective function 2 of 



Sensors 2009, 9 



8696 



problem P3 aims at finding the best location for sp.^-^that minimizes the transmission loss of the link 
l(i-x)k-+\(i-x+i)k-+\ or the link l(i- x -\)k (i-x)k-^i with the highest transmission loss value and preserves the 
connectivity with node S(/- x +i)^i and nodes £(/_ x _i)£ preventing from network partition. Objective 1 is 
valid when the nodes S(i- x -i)k are also located in the shadow zone. Otherwise objective 2 should be 
applied. Constraints (18) and (19) impose that the transmission range should be lower than a threshold 
to ensure connectivity. Constraints (20) and (21) express that the transmission loss should be lower 
than a transmission loss threshold. Constraints (19) and (21) are only valid for objective function 2. 
Constraint (22) imposes that the sensor node S(i- x )k-+i should be located outside the shadow zone. 
Constraint (23) states that the energy consumed by the electronically controlled engine to move the 
sensor to the optimal position should be lower than a threshold. 



Now we study the performance of the proposed scheme under the presence of shadow zones via 
numerical evaluations. We distinguish between shallow and deep water. The five different channel 
models considered appear in Section 4. The parameters used in our evaluation are listed in Table 1. 
They follow the architecture of the commercial off-the-shelf underwater acoustic modem (LinkQuest) 
UWM3000 for shallow and UWM10000 for deep water [22]. 



5. Results 



Table 1. Parameter Values. 



Parameter 



Value 



Sea depth H 



Shallow water: 100 m 
Deep water: 5,000 m 



Volume 



Shallow water: 16,000 x 16,000 x 100 m 
Deep water: 

40,000 x 40,000 x 5,000 m 3 



Number of sensors 
Packet size 
Data rate 
T 



25 

238 bytes 
6 Kbps 
15 °C 



pH 
S 



35 ppt 

Shallow water: 12 W 
Deep water: 40 W 
47.69 dB 
1500 m/s 
0.5 m/s 



8 



Transmission power P t 



Noise level 
Speed of sound 
Sensor speed 
Movement power P Bw 

A 



0.01 



10W 



Rmax 



TL Th 



Shallow water: 5,000 m 
Deep water: 1 1,000 m 
1,500 J 

Shallow water: 112.5 dB 
Deep water: 129 dB 
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The whole seabed is divided into a 2D square grid of equal sizes r = 4,000m for shallow water and 
r = 10,000 m for deep water as shown in Figure 6. The sensor nodes are located in a square grid 
floating at different depths. 25 sensor nodes are deployed in a 3D volume of 16,000 x 16,000 x 100 m 3 
for shallow water and of 40,000 x 40,000 x 5,000 m for deep water. In shallow water, sensor nodes 
are located following a hierarchical structure at the depth levels of 25, 50, 75 and 90 m. In deep water, 
we study the communication between sensor nodes located following a hierarchical structure at the 
depth levels of 300, 350, 400 and 450 m. 



Figure 6. 2D Square grid. 

y 




2DXY- 
Square grid 




Sensor nodes placed in the 



r ■ 10000 m (deep water) square grid floating at different 

depths 

We consider realistic cases where shadow zones can influence the propagation of sensors. 
Figures 7, 8 and 9 show examples of the propagation paths for the sensor nodes in shallow water. The 
sound speed profile is based on the measurements conducted at the coast of Kauai, Hawaii [6], with a 
water depth of 100 m and a frequency of 1 kHz. The rays shown are for launch angles between -20° 
and 20° with a 1° increment. The corresponding ray traces in Figure 7 for a source depth of 90 m and a 
receiver depth of 75 m, in Figure 8 for a source depth of 75 m and a receiver depth of 50 m and in 
Figure 9 for a source depth of 50 m and a receiver depth of 25 m are plotted. Shadow zones appear at 
ranges around 3 Kms. In our simulations we consider a 2D square grid of equal sizes r = 4,000 m 
because at these ranges it is possible to appreciate the shadow zone effects. 
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Figure 7. Corresponding ray trace for a source depth of 90 m. 



BELLHOP- f=1000Hz 




Range (m) 

Figure 8. Corresponding ray trace for a source depth of 75 m. 



BELLHOP- f=1000Hz 




Range (m) 



Figure 9. Corresponding ray trace for a source depth of 50 m. 
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Figures 10, 11 and 12 show examples of the propagation paths for the sensor nodes in deep water. 
We consider a frequency of 1 kHz and a deep-water (5,000 m) environment with a Munk sound speed 
profile [2]. The rays shown are for launch angles between -20° and 20° with a 1° increment. The 
corresponding ray trace in Figure 10 for a source depth of 450 m and a receiver depth of 400 m, in 
Figure 11 for a source depth of 400 m and a receiver depth of 350 m and in Figure 12 for a source 
depth of 350 m and a receiver depth of 300 m are plotted. Shadow zones appear at ranges around 10 Kms. 
In our simulations we consider a 2D square grid of equal sizes r = 10,000 m because at these ranges it 
is possible to appreciate the shadow zone effects. 

Figure 10. Corresponding ray trace for a source depth of 450 m. 

BELLHOP- f=1000Hz 




Figure 11. Corresponding ray trace for a source depth of 400 m. 

BELLHOP- f=1000Hz 
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Figure 12. Corresponding ray trace for a source depth of 350 m. 

BELLHOP- f=1000Hz 




Transmisson Loss Threshold 

Now we compute the transmission loss threshold TL T h or maximum propagation loss for properly 
receiving the transmitted signal given by (9). The link quality in underwater communication is 
severely affected by multipath phenomena, multipath fading and the refractive properties of the sound 
channel [3]. As a result, the bit error rates of the acoustic links are often high and efficient error control 
schemes are needed. The error control schemes automatic repeat request (ARQ), forward error 
correction (FEC) block codes and FEC convolutional codes have been selected and evaluated in terms 
of the packet error rate (PER). 

According to [3], OFDM transmission with QAM modulation has been used for our analysis. The 
scalable OFDM with 16-QAM modulation has been successfully tested in the settings with 
bandwidths 12 kHz, 25 kHz and 50 kHz, leading to data rates about 12 kbps, 25 kbps and 50 kbps 
respectively [23]. The BER for the modulation scheme 16-QAM is given by: 



Pb 



16QAM _ 



8 



erfc 



4 E u 



10N, 



E b IN 0 = ¥ 



o J 



R 



(24) 



where ^lO^ 5 ^ 710 , B N is the noise bandwidth and R is the data rate. 

Based on the bit error rate pi» the PER for the error control schemes can be calculated as follows: 
For ARQ, the cyclic redundandy check (CRC) block code detection mechanism is deployed. 

Assuming detection of all possible packet errors, the PER of a single transmission for a packet of / bits 

is computed as: 

PER CRC (l) = l-{l-p b ) 



(25) 



A BCH (Bose, Ray-Chaudhuri, Hocquenghem) and a RS (Reed-Solomon) FEC block code is 
represented by (n, k, t), where n is the block length, k is the payload length, and t is the error correcting 
capability in bits. For the BCH and RS codes, the block error rate (BLER) is given by: 



BLER(n,k,t)=Z . P%~P b V 



i=t + l 



(26) 



Since a packet can be larger than the block length n, the PER for BCH and RS block codes is given by: 
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PER BC (l,n,k,t) = l-(l-BLER(n,k,t)) * (27) 

where \l lk \ is the number of blocks required to send / bits and is the ceiling function. 

For FEC convolutional codes, the PER of a single transmission for a packet of / bits is given by: 

\i] 

pER convolutional ^ Wj jfc) = 1 - (l " /?J ° (28 ) 

where pb is the bit error probability of an encoded data packet of length of length \l /R c ~] and code 
rate R c =k/n. 

The relationship between the PER and TL T h for ARQ, three BCH, three RS and two convolutional 
codes is shown in Figure 13 for shallow water and in Figure 14 for deep water. For a target PER, FEC 
block codes can support TL T h values higher than ARQ and convolutional codes. We notice that the 
maximum TL T h value is increased with the error correcting capability of the FEC block codes and the 
code rate of the FEC convolutional codes. For a particular TL T h value, the PER is equal or lower for 
error control codes in deep water in comparison with shallow water. 

Figure 13. PER as a function of transmission loss for shallow water. 
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Figure 14. PER as a function of transmission loss for deep water. 
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The FEC convolutional code with R c = 1/2 has the lowest TL T h value. This is the worst case, 
because the maximum transmission loss is restricted to lower values. Therefore, with this error control 
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scheme for a target PER of 10 the maximum TL T h of 1 12.5 dB for shallow water and 129 dB for deep 
water have been selected. More details about the performance of error control schemes in UWSNs can 
be found in [20]. 

Energy Consumption in the Sensor Movement 

Now we analyze the energy consumed by the electronically controlled engine to move the sensors 
affected by the shadow zone to the optimal position. Figure 15 shows the energy consumption in the 
sensor movement. We observe that the energy consumed is increased with the movement power and 
the covered distance to reach the optimal location. Since the sensor relocation can be costly in terms of 
energy consumption, it is limited by the threshold E T h\ this paramenter is estimated analyzing the 
battery capacity, the required average displacement as well as P Bw - 

Figure 15. Energy consumption in the sensor displacement. 

^ _ ^ - - - ' r ' ' 
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Case 1: Optimal Placement Evaluation 

We consider a sensor node s t located at y i9 z t ). Node s i+ \ is located at (jc,-± j.r, yt± j.r, z i+ i),j e N, 
r = 4,000 m for shallow water, r = 10,000 m for deep water, z i+ \ < z t . The acoustic link lu+\ is affected 
by a shadow zone with height (z sz - end ~ z sz - start ) = 9 m, where z t > z sz - end and z i+ \ < z sz - start . Node s t is 
uncoupled into sensor nodes s^i and s z -_>2 and the optimal location of node s z -_>2 has been found using 
NLP solving the optimization problem PI. 

Figure 16 shows the average transmission loss as a function of frequency for shallow and deep 
water under different propagation phenomena. We have computed the transmission loss for the 
acoustic link l ii+ i (between s t and s i+ \) affected by the shadow zone (shallow water + shadow zone or 
deep water + shadow zone). St is uncoupled and the optimal location for node s^2 is found outside the 
shadow zone using our mathematical model; the transmission loss for the acoustic link /;_> 2 i+i 
(between 5,-_>2 and s i+ \) (shallow water (optimal), deep water + convergence zone (optimal) or deep 
water + deep sound channel (optimal) has also been computed. The average transmission loss is 
increased when the frequency is increased and is higher for deep water + shadow zone than for shallow 
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water + shadow zone. It is also higher for deep than shallow water. The transmission loss values are 
higher than the transmission loss threshold for deep water+shadow zone, which means that the 
transmitted signal will not be properly received. The average transmission loss values are reduced 
significantly when the optimal locations are computed and the communication between sensor nodes is 
again restablished outside the shadow zone. The maximum transmission loss improvement is of 4.9 dB 
for shallow water, of 45.1 dB for deep water + convergence zone and of 42.9 dB for deep water + deep 
sound channel. Therefore, we can conclude that using our mathematical model deep water + convergence 
zone shows the best improvement in the diminishment of the transmission loss. 

Figure 16. Average transmission loss of a link as a function of frequency in shallow and 
deep water. 




0 2 4 6 8 10 12 

Frequency (kHz) 



Case 2: Optimal Placement Evaluation 

We consider a sensor node s\ located at (jc,-, yu zi). Node s i+ \ is located at (jc,-± j.r, y x ? ± j.r, z i+ i),j e N, 
r = 4,000 m for shallow water, r = 10,000 m for deep water, z i+ \ < z t . Nodes S(i-\)k, V£ e K are located 
at (xt±j.r, yt±j.r, Z( z _i)^),y e N, r = 4,000 m for shallow water, r = 10,000 m for deep water, z^k > z t , 
k g K. The acoustic links lu+i and l(i-\)ki are affected by a shadow zone with height (z sz -end ~ Zsz-stan), 
where z sz - start < z t < z sz - end , z (iA)k > z sz - end , z m < z sz - start . Node s t is uncoupled into sensor nodes s^i and 
5,-_>2. The optimal locations of node s^2 and ^i have been found using NLP solving the optimization 
problems PI and P2, respectively. Now we evaluate the performance of the proposed scheme for the 
optimal location of node s^\ (Case 2). For shallow water the value of z sz - en d varies between 50 and 85 
meters, z sz - start = 40 m andz^.i)* e (z^-^lOOJ.The average transmission loss of the acoustic link l(i-\)ki 
as a function of the shadow zone height (z sz - en d - z sz - sta rt) for shallow water and three different 
frequencies has been computed. The transmission loss of the link l(i-\)ki between the node s^k, Vk e K 
and Si^i with the highest value is minimized for shallow water. The transmission loss values are not 
affected by the shadow zone height but they are increased with the frequency to 41 dB for 
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/=0.1 kHz, 41.6 dB for f=\ kHz and 49.3 dB for /=10 kHz. Figure 17 shows the average transmission 
loss of the acoustic link l(i.\)ki as a function of the shadow zone height (z sz - en d - z sz - start ) for deep water 
and three different frequencies. For deep water the value of z sz -end varies between 500 
and 4000 meters, z sz - start = 440m and z^k e (z^-^,5,000]. The transmission loss of the link ,■ 
between the node S(i-\)h Vk e K and s^\ with the highest value is minimized for different propagation 
phenomena (deep water + deep sound channel, deep water + convergence zone). Deep water + 
convergence zone is the underwater propagation phenomena that suffers lower transmission loss and is 
more appropriate for underwater communication. 

Figure 17. Average transmission loss of a link as a function of the shadow zone height in 
deep water. 
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The transmission loss values are especially decreased with the shadow zone height for the 
frequency of 10 KHz; the reason is that the depth of the optimal location is increased when the shadow 
zone height is increased (z sz - en d is higher) and the transmission loss decreases when the depth is 
increased. Otherwise, the transmission loss values for 0.1 KHz and 1 KHz are only very slightly 
decreased by the shadow zone height and are only significantly increased with the frequency. 

Case 3: Optimal Placement Evaluation 



We consider a sensor node Si located at (jc,-, yi, Zi). Nodes S(/_i)yb Vk e K are located at (x z ± j.r, 
yi± j.r, Z(i-i)k) 9 j e N, r = 4,000 m for shallow water, r = 10,000 m for deep water, , z^iyt > z u k e K. 
Nodes S(i-2)k, \/k e K are located at (jc,-± j.r, yt± j.r, Z( z _2)^),y e N, r = 4,000 m for shallow water, r = 
10,000 m for deep water, , z^ 2 )k > z (i-i)h k e K. Nodes s^yh Vk e K are located at (x t ± j.r, y t ± j.r, 
Z(i-3)k),j g N, r = 4,000 m for shallow water, r = 10,000 m for deep water, z^t > z^k, k e K. The 
acoustic links l^.\yku h-iykii-vyk , l(i-3)k(i-iyk are affected by a shadow zone with height (z sz - end - z sz - start ), 
where z sz - sta rt < z\ < z sz - en d, z sz - sta rt < ^(z-i)A: < z sz - en d, z sz-start < z (i-2)k < z sz - en d, Z(j-3)k > z sz - en d. Node Si is 
uncoupled into sensor nodes s^i and s^2- Nodes s^iyk, \/k e K are uncoupled into sensor nodes 
S(i-iyt-+i and S(i-iyt-+2- Nodes s^k, Vk e K are uncoupled into sensor nodes s^k^i and s^k^i- Node 
s^\ is relocated to (x z , y t , z sz _ en d+s), where s e i? + is a very small value, z sz - en d= 89 m for shallow water 
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and z sz - en d = 4049 m for deep water. The optimal locations of nodes S(i-\yt-+i have been found with NLP 
solving the optimization problem P3 with the objective function 1. The optimal locations of nodes 
S(i-2)k-+i have been found with NLP solving the optimization problem P3 with the objective function 2. 
In shallow water, sensor nodes are located following a hierarchical structure at the depth levels 
of 25, 50, 75 and 90 meters. In deep water, sensor nodes are located following a hierarchical structure 
at the depth levels of 3,900, 3,950, 4,000 and 4,050 meters. 

We consider four hierarchical levels with one node Si at the lowest level, three nodes SQ-iyt, Vk e K 
at the next level, six nodes s^yk, V£ e K at the next level and 12 nodes s^yt, V£ e K at the deepest 
level. Figure 18 shows the average SNR for the acoustic links l(i-\)ku ki-i)k(i-\)k, l(i-3)k (i-i)k affected by a 
shadow zone as a function of frequency for shallow and deep water under different propagation 
phenomena. The average SNR values in the presence of shadow zones are negative, which indicates 
that the signal is below the noise level and transmission loss, and it can't be properly recovered. They 
have not been depicted but are decreased when the frequency is increased and are lower for deep 
water + shadow zone than for shallow water + shadow zone because the transmission loss values are 
higher for deep water + shadow zone. The average SNR values are increased significantly when the 
optimal locations are computed and the communication between sensor nodes is again restablished 
outside the shadow zone. Furthermore, the improved SNR values are maintained for all frequencies. 
With shallow water the SNR values are better because in shallow water the transmission loss values 
are lower. The SNR values are higher with deep water + convergence zone than with deep water + deep 
sound channel. Nevertheless, in all cases the SNR values are positive and they are increased outside 
the shadow zones. The maximum SNR improvement is very significant: 197.64 dB for shallow 
water, 352.08 dB for deep water + convergence zone and 350.68 dB for deep water + deep 
sound channel. 

Figure 18. Average SNR as a function of frequency in shallow and deep water. 
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6. Conclusions 

In this paper, a distributed adaptive topology reorganization scheme that maintains connectivity in 
multi-hop UWSNs affected by shadow zones has been developed in the context of two Spanish-funded 
research projects. It alleviates the effect of energy limitations and is able to maintain network 
connectivity in multi-hop three-dimensional UWSNs in the presence of shadow zones solving a 
location optimization problem. 

Three different cases have been determined according to the number of links in the three-dimensional 
UWSN affected by the shadow zone. A mathematical model has been developed for each case to find 
the optimal placement for the sensor nodes, whose ongoing communications are being disturbed. 

The average transmission loss values are reduced significantly (especially for deep water + convergence 
zone) when the optimal locations are computed and the communication between sensor nodes is again 
reestablished outside the shadow zone. Shallow water shows lower transmission loss values in 
comparison with the other propagation phenomena. Deep water + convergence zone suffers lower 
transmission loss values than deep water + deep sound channel and is more appropriate for underwater 
communication. Only in deep water the transmission loss values are especially decreased with the 
shadow zone height for the frequency of 10 KHz and only very slightly decreased for the frequencies 
of 0.1 KHz and 1 KHz; the transmission loss values don't vary with the shadow zone height in 
shallow water. 

The average SNR values in the presence of shadow zones are negative, which indicates that the 
signal is below the noise level and transmission loss, and it can't be properly recovered. The average 
SNR values are increased significantly when the optimal locations are computed and the 
communication between underwater sensors is again reestablished outside the shadow zone. 
Furthermore, the improved SNR values are maintained for all frequencies. With shallow water the 
SNR values are better because in shallow water the transmission loss values are lower. The SNR 
values are higher with deep water + convergence zone than with deep water + deep sound channel. 
Nevertheless, in all cases the SNR loss values are increased outside the shadow zones. 
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